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ABSTRACT 

Results are presented from a survey for molecular anions in seven nearby Galactic star-forming cores and 
molecular clouds. The hydrocarbon anion C<,H" is detected in all seven target sources, including four sources 
where no anions have been previously detected: LI 172, L1389, L1495B and TMC-1C. The C6H"/CeH column 
density ratio is > 1 .0% in every source, with a mean value of 3.0% (and standard deviation 0.92%). Combined 
with previous detections, our results show that anions are ubiquitous in dense clouds wherever C(,H is present. 
The CeH"/C6H ratio is found to show a positive correlation with molecular hydrogen number density, and 
with the apparent age of the cloud. We also report the first detection of C4H" in TMC-1 (at 4.8er confidence), 
and derive an anion-to-neutral ratio C 4 H7C 4 H = (1 .2 ± 0.4) x 10" 5 (= 0.0012 ± 0.0004%). Such a low value 
compared with CeH~ highlights the need for a revised radiative electron attachment rate for C4H. Chemical 
model calculations show that the observed C4H~ could be produced as a result of reactions of oxygen atoms 
with C 5 H" and C 6 H". 

Subject headings: astrochemistry — ISM: abundances — ISM: molecules — ISM: clouds — stars: formation 



1. INTRODUCTION 

Negative ions (anio ns) were first discovered in the interstel- 
lar medium (ISM) by iMcCarthy et al.1 (120061) . who detected 
the linear hydrocarbon anion CeH" in TMC-1. This was fol- 
lowed by CeH~ and C^H" de tecti ons in the protostellar core 
L1527 by ISakai et al.l d2007l) and lAgundez et aTJ (120081) . re- 
spectively, which hint ed of a role for anions in the chem- 
istry of star-formation. Gupta et al. (2009) performed the first 
dedicated survey for CeH" in twenty-four molecular sources, 
but detected the anion in only two star-forming clouds, with 
anion-to-neutral ratios on the order of a few percent, similar 
to previously-observed values. Interstellar anion detections 
were confined to the Taurus Molecular Cloud complex un- 
til the recent discoveries of C6H~ in the Lupus, Cepheus and 
Auriga star-forming regions (Sakai et al .120101 ICor diner et al. 
1201 11) . which proved that anions are widespread in the local 
ISM, and not an artifact of any particular physical or chemical 
conditions in the Taurus region. 

The possible importanc e of anions in interstellar chem- 
istry was first discussed by Dalgarno & Mc Cravl (11973b . and 
prospects for the detection of molecular an io ns usin g radio 
astronomy were examined by iSarrel (11980b . iHerbsj (11981b 
argued that large interstellar molecules (including carbon- 
chain-bearing species) can undergo rapid radiative electron 
att achment - as exem plified by the laboratory experiments 
of iWoodin et al.l (119801) - potentially resulting in significant 
anion abundances in dense molecular clouds. Based on this 
idea, models for anion chemistry have been successful in re- 
producing the observed abundances of C(,H~ and CgH" in 
TMC-1, IRC+10216, L1527 and L1512 Millar et d J 120071: 
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Remii an et ail 120071: lHarada & Her bst 2008; ICordiner et aTJ 
2008: ICordiner & Charnlevll2012b . However, there is a ma- 
jor discrepancy between the modeled and observ ed C4H" 
anion-to-neutral ratio (see Herbst & Osamura 2008), and the 
lack of anions in PDRs ( Agundez et al. 2008) is at variance 
with the model predictions of Millar et al. (2007). Clearly, 
our understanding of molecular anion chemistry is incom- 
plete. Nevertheless, as demonstrated by ICordiner & Charnlevl 
d2012l) . anion measurements have the potential to offer in- 
sight into interstellar clou d properties due to their reactivity 
( Eichelber ger et al.l 120071) . and consequently, their sensitiv- 
ity to the abundances of gas-phase electrons and atomic C, 
H and O, and to physical conditions such as cloud density 
and molecular depletion. Observations of anions in a variety 
of interstellar environments will be key to a complete under- 
standing of their chemistry. Given the relatively small number 
of molecular anion detections in the ISM to date, we set out 
to address the question of just how widespread anions are in 
interstellar clouds, and to ascertain the behavior of the anion- 
to-neutral ratio over a range of cloud types and ages. 

In this article, results are presented from our survey for 
CgH" and CeH in a sample of seven carbon-chain-rich inter- 
stellar clouds, protostars and prestellar cores. Detections of 
CgH" in two of th e sources (LI 251 A an d L1512) were previ- 
ously reported by Cordi ner" et al.l (1201 lb . and here we report 
the final results for all seven sources, as well as the first (ten- 
tative) detection of C 4 H" in TMC-1. 

2. TARGET SOURCES 

Due to the close chemical relatio nship between p o lyynes 
and cyanopolyynes (see e.g. Fede rman et al.l 119901: 
Millar & Herbst 1 19941) . carbon-chain-bearing species such as 
C4H and CeH are expected to be abundant in dense molecular 
clouds, close to the peak of HC3N emission. Target sources 
with strong HC3N emission lines were selected from an HC3N 
J = 10-9 mapping survey of twenty nearby molecular clouds, 
prestellar cores and young, low-mass protostars obtained us- 
ing the Onsala Space Observatory (OSO) 20-m telescope be- 
tween 2005 and 20 12. S ome of these data wer e presented by 
iBuckle et all d2006l) and ICordiner et al l (12011b . Maps of the 
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TABLE 1 

Targets, coordinates and distances 



Source 


RA 


Dec 


Cloud Type" 


Distance 


Ref. 




(J2000) 


(J2000) 




(PC) 




LI 172 SMM 


21:02:22.1 


+67:54:48 


Star-forming 


440 


1 


L1251A 


22:30:40.4 


+75:13:46 


Star-forming 


300 


2 


L1389 (CB17) SMM1 


04:04:36.6 


+56:56:00 


Protostellar 


250 


3 


L1495B 


04:15:41.8 


+28:47:46 


Quiescent 


140 


4 


L1512 


05:04:07.1 


+32:43:09 


Star-forming 


140 


4 


TMC-1C 


04:41:35.6 


+26:00:21 


Star-forming 


140 


4 


TMC-1 CP 


04:41:41.9 


+25:41:27 


Quiescent 


140 


4 



REFERENCES. — (1) Visser et al. 12002); (2) Kun & Prusti 1 1993); (3) Launhardt et al. (2010); (4) Elias ( 1978). 
a See Section |2]for description of source classification scheme. 



integrated HC3N J = 10-9 emission intensity towards six of 
our seven chosen sources (LI 172, L1251A, L1389, L1495B, 
L1512 and TMC-1 C) are shown in Figure [TJ These observa- 
tions were obtained using the OSO with a beam size of 42" 
and beam efficiency of 0.5 ±0.05. A 30" map spacing was 
used, which was resampled to a 7.5" per pixel grid using bi- 
lin ear interpolation. For TM C-1, HC3N maps we re published 
bv lHirahara etaTI (fl992b and lPratap et afl(U997l) . 

To search for CeH and CeH", we targeted the strongest 
HC 3 N peak positions in our maps of L1389, L1495B, L1512 
and TMC- 1C. For LI 172 we targeted the sub-mm core cat- 
aloged by iDi Francesco et al.l d2008l) . which fell within 15" 
of the HC3N peak. Our chosen LI 251 A position does not 
coincide with the main HC3N peak because it was based on 
an earlier, lower-sensitivity map than that shown in Figure Q] 
For TMC- 1 , w e targe ted the cyanopolyyne peak (CP) where 
Agun dez et al.1 (120081) previously attempted to detect C4H". 
The adopted coordinates for our anion survey targets are listed 
in Table □ 

Our LI 172 and L1512 Ce H" positions overla p the 
prestellar cores iden ti fied by IVisser et"ai1 (120021) and 
IWard-Thompson et al.l (119941) . respectively. For L1389, 
the HC3N emission maps out a very compact clump, 
coincident with the cente r of the Bok Globule CB17 
(IClemens & Barvainislll988l) . The HC3N peak also matches 
closely the location of strong est microwave/ sub-mm contin- 
uum emission observed by Laun hardt et al.l (120101) and de- 
noted L1389 SMM1 (the position of which is indicated in the 
bottom-left panel of Figure [TJ by the cross inside the dashed 
black circle of the GBT beam). As shown in Figure [TJ our 
targeted L1251A position is ~ 40" south- east of the center of 
the protostar L125 1 IRS3 dLee et al.l2010l) . although the outer 
envelope of this protostar prob ably intersects the GB T beam 
to some extent, as discussed bv lCordiner et al.l (1201 ll). 

From dust continuum emission, Launha rdt et al J fob 101) ob- 
tained a number den s ity for L1389 SMM1 of 6 x 10 6 cm" 3 . 
iPavlyuchenkov etaf] (120061) derived a moderately high CO 
depletion factor (~ 40), and inferred a relatively evolved 
chemical age of ~ 2 Myr for this collapsing core. The core 
was recently found to contain a very young, low-lumin osity 
protostar (referred to as CB17 MMS; IChen et al.l [20121) . 8" 
NE of our targeted GBT beam center, as well a more evolved 
Class 0/1 protostar 20" NW. Given that the radii of typ- 
ical low-mass protos tellar envelopes are ~ 10 4 AU {e.g. 
Jorge nsen et al.ll2002l) . at a distance of 250 pc, our targeted 
L1389 beam is likely to be dominated by protostellar enve- 
lope matter. 

L1495B, on the other hand, is a quiescent molecular cloud 
with no known protostars nearby. It appears to be chemi- 



cally young, with an apparently low level of C O depletion 
consi stent with less- evolved inters t ellar c louds dHirota et alj 
I2004T) . In TMC-1C. Iluckle et alj d2006h identified an anti- 
correlation between the spatial distributions of HC3N and 
C ls O, at least partly attributable to CO depletion. Thus, our 
observations of TMC-1C (at the HC3N peak) probably sam- 
ple more chemically-evolved, depleted gas than L1495B. No 
protostars are known to be present in our observed TMC- 
1C beam, but several sub-mm sources are located in the sur- 
rounding cloud (indicated by white crosses in the lower-right 
panel of Figure [T), showing this to be a region of active star 
formation. South-west of TMC-1C lies the source TMC-1 
(CP), which is a well-known chemical ly-young dark cloud , 
with only a modest degree of depletion (Hirahar a et al.lll 992t 
ICordiner & Charnlevll2012l) . 

In Table [TJ we provide a basic categorization of our target 
sources in light of their properties described above: 'quies- 
cent' refers to those cloud cores that are chemically young, 
show no evidence for active star-formation within the tele- 
scope beam (such as outflows or compact sub-mm/IR emis- 
sion), and do not appear to be undergoing collapse; 'star- 
forming' is used for clouds showing nearby active star- 
formation and a greater degree of chemical evolution; finally, 
'protostellar' describes our L1389 position, which contains 
the low-luminosity protostar CB17 MMS. 

3. ANION OBSERVATIONS 

Observations of emission lines of CeH (J = 10.5-9.5), 
C 6 H" (J = 10-9 and 11 - 10) and HC 3 N (J = 3-2 and 4-3) 
were carried out between 2010 April and 2011 January us- 
ing the Ka receiver of the NRAO 100-meter Robert C. Byrd 
Green Bank Telescope^ (GBT). The GBT Spectrometer was 
used with a bandwidth of 50 MHz and 8192 x 6.1 kHz chan- 
nels (corresponding to a velocity spacing of « 0.065 km s" 1 ). 
Four spectral windows were used, which allowed simultane- 
ous observation of the CeH and CeH" lines. For the compact, 
spatially-isolated source L1512, beam switching (with 78" 
throw) was used, and for all other targets, frequency switch- 
ing was used. Pointing was checked every one to two hours 
and was typically accurate to within 5". In the middle of 
the observed frequency range (28 GHz), the telescope beam 
FWHM was 26" and the main beam efficiency was 0.90. To- 
tal system temperatures were typically in the range 60 - 80 K 
and the zenith opacity was 0.05 ± 0.02. Intensity calibration 
was performed using beam-switched observations of the com- 
pact radio source NGC 7027. Measured antenna temperatures 

5 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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FIG. 1.— Observed HC 3 N J = 10-9 integrated intensity maps towards L1172, L1251A and L1389, L1495B, L1512 and TMC-1C. The OSO beam size is 
indicated in the lower right of each panel with a white circle. Intensity scales are given in units of Kkms~' and have not been corrected for beam efficiency. 
Dashed black circles represent the size (HPBW) and position of our targeted GB T anion survey po sitions. Crosses denote locations of sub-mm (prestell ar) cores 
jWard-Thompson et alj fl994; Visser et al. 2002: T Di France ' ' " 

2002; Chen et al. 2012; Lee et al. 2010). 



i Francesco et al. 2008; Launhardt et al. 2010) and asterisks denote protostellar core locations jVisser etaT 
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were subsequently corrected for opacity, spillover, ohmic loss, 
blockage efficiency and beam efficiency, then averaged using 
standard GBTIDL routines. 

We obtained observations of C 4 H (N = 2- 1, J = 1.5-0.5, 
F = 2 - 1) and C4H" (7 = 2-1) between 201 1 December and 
2012 February at around 19 GHz using the GBT K-band 
focal-plane array (KFPA). These observations were obtained 
with channel spacing 12.2 kHz, beam efficiency 0.92, beam 
size 39" and zenith opacity 0.025 ± 0.015. 

Additional single-pointing HC3N J = 10-9 spectra of our 
sources were obtained using the OSO at our adopted target 
positions. A list of all the observed transitions is given in 
Table |2] 

4. RESULTS 

4. 1 . New Anion Detections 

The CeH" anion and its neutral counterpart CeH were de- 
tected in all of the surveyed sources. This is the first time an- 
ions have been detected in LI 172, L1389, L1495B and TMC- 
1C. The observed CeH" and CeH spectra are shown in Figure 

m 

Spectra of the J = 10-9 and J = 11-10 transitions of C^W 
were averaged in velocity space to improve the signal-to-noise 
ratio; integrated intensities of the averaged spectra are given 
in Table [3] Due to partial blending of the hyperfine compo- 
nents of the CeH (J = 10.5-9.5,/) lines, the integrated inten- 
sity summed over both components is given in Table [3] Given 
the weakness of these lines, they ca n safely be ass umed to be 
optically thin, and Equation 2 of iLis et all (2002) was used 
to derive column densities for CeH and CeH" (given in Table 
0. Spectroscopic data for the transitions of interest (Tabled 
were obtai ned from the Colo gne Database for Molecular As- 
tronomy (Miille r et al.l 120051) . and an excitation temperature 
of 10 K was assumed. Error estimates were derived using 500 
Monte Carlo noise replications for each measurement. 

In TMC-1, CeH excitation temper ature s of 5.2 K and 6.7 K 
were measured by iBell et al.1 (11999b and ISakai et al.l (120071) . 
respectively. Such sub-thermal excitation might be expected 
given the large Einstein A coefficients for rotational transi- 
tions of this molecule, which are a direct consequence of its 
lar ge (5.5 D) dipole moment. However, the value obtained 
bv Bell et al.l (Tl999) may be subject to un certainty due to the 
effects of telescope beam dilution, and the Sakai et al. (2007) 
value may not be applicable to our observations because it 
was obtained using the J = 16.5- 15.5 transition, which has 
an Einstein A coefficient about four times greater than that of 
the /= 10.5-9.5 transition we observed. The larger Einstein 
A results in an increased likelihood of sub-thermal excitation 
of the J = 16.5 level. If the CeH excitation temperature is as 
low as 5 K, then the calculated column densities for CeH and 
CeH" will be about 36% smaller than the values in Table [3] 
Due to the similar moments of inertia of CeH and CeH", their 
rotational levels occur at similar energies, so that the calcu- 
lated anion-to-neutral column density ratios (CeH"/CeH) are 
insensitive to small uncertainties in excitation, provided both 
molecules share a common excitation temperature. The mean 
CeH"/CeH ratio for the seven sources is 3.0%, with a standard 
deviation of 0.92%. 

Our new measurements for CeH an d CeH" in L125 1A an d 
L1512 supersede those presented by Cordi ner et al.l (1201 lb . 
Although the column densities and anion-t o-neutral ratios 
measu red in the present study match those of ICordiner et al.l 
d201 ll) within the stated errors, the newer values benefit from 



improved calibration and analysis methods and have smaller 
uncertainties. The difference between our TMC-1 (CP) 
CeH-/C 6 H value of 2.5 ± 4% a nd the value of 1.6 ±0.3% 
obtained by Briinke n'et al.l (120071) is just beyond the range of 
the one-sigma error bars, and could be the re sult of statisti- 
cal no ise in the spectra. Alternatively, because Briinken et al. 
(120071) derived their value from a different set of transitions, 
non-LTE excitation effects in CeH and/or CeH" could be par- 
tially responsible for the discrepancy. 

The observed C 4 H" and C 4 H spectra of TMC-1 (CP) are 
shown in the bottom-right panel of Figure [2] The 7 = 2-1 
line of C4H" is detected in a single channel with a peak 
antenna temperature of 3.4 mK, which corresponds to 4.8<r 
(where a is the RMS noise of the baseline-subtracted spec- 
trum). Using a Gaussian fit to this C 4 H~ line with 1000 
Monte Carlo noise replications, the central velocity was 
found to be 5.70 ±0.05 kms" 1 with a FWHM of 0.43 ± 
0.13 kms" 1 . Within the errors, these parameters match those 
derived for HC3N (given in Table |4), which adds confi- 
dence to our C 4 H" detection. The integrated line intensity 
is 1.0 ±0.3 mKkms" 1 , which corresponds to a column den- 
sity of Af(C 4 H~) = (8.0 ± 2.4) x 10 9 cm" 2 (assuming an exci- 
tation temperature of 10 K). The neutral C 4 H line has an in- 
tegrated intensity of 41 1.3 ±0.7 mKkms" 1 , which requires 
JV(C 4 H) = (6.91 ±0.01) x 10 14 cm" 2 . The corresponding 
anion-to-neutral ratio C 4 H"/C 4 H is (1.2 ± 0.4) x 10" 5 . This 
value is consistent with the u pper limit of 5.2 x 10~ 5 mea- 
sured by Agundez et al. (2008), and our C 4 H column density 
also closely matches their value of 7.1 x 10 14 cm" 2 . Both 
Af(C 4 H) and A^(C 4 H~) are relatively insensitive to the adopted 
excitation temperature, and vary by only 12% over the range 
5-10K. 

4.2. HCjN Observations and Radiative Transfer Modeling 

Our HC3N maps (Figure [TJ show the presence of one or 
more compact molecular condensations in each source, and 
highlight a clumpy structure in this carbon-chain-rich gas. 
The peaks in HC3N emission tend to coincide approximately 
with the locations of known prestellar cores and protostars, 
the positions of which are denoted with crosses and asterisks, 
respectively. As a result of the collapse and infall processes 
occurring in these objects, their densities are expected to be 
greater than the surrounding gas. This should result in in- 
creased collisional excitation of the J = 10 le vel of HC3N, 
which has a critical density of 7 x 10 5 cm" 3 dBuckle et alj 
2006). Thus, the intensity variations shown in ourHC3N OSO 
maps represent a combination of variations in gas density and 
total HC3N column. 

Spectra of three rotational transitions of HC3N for each 
source, observed at the locations targeted by our GBT an- 
ion survey, are shown in Figure [3] These span a range of 
upper-state energy levels 1 .82- 16.69 cm" 1 (2.6-24.0 K). In 
order to derive gas densities and velocities, the HC3N spec- 
tra were subject to fitting using the RADEX radiative trans- 
fer code (Ivan der Tak et al.ll2007l) . Collisional transition rates 
were taken from the Leiden Atomic and Molecular Database^ 
(Schoi eFet al.ll2005|). which tabulates sca led versions of the 
original data of iGreen & Chapmanl (119781) . Hyperfine struc- 
ture (HFS) was accounted for by assuming LTE population 
of the hyperfine levels in each J state. The primary collision 
partner number density («h,), molecular column density (AO, 

6 http://www.strw.leidenuniv.nl/~moldata 
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FIG. 2. — Observed CgFT spectra (average of J = 10 — 9 and / = 11-10 transitions), and CgH / = 10.5 — 9.5, / spectra. The CgH velocity scale is given with 
respect to the (weighted) mean frequency of the two hyperfine components. For TMC-1, observed Cz|I-r {J = 2—1) and C4H (N = 2—1,7 = 1.5 — 0.5, F = 2—1) 
spectra are also shown (bottom-right panel). Dashed vertical lines indicate the least-squares HC3N velocities for each target. 
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TABLE 2 

Observed species, transitions and frequencies 



Species 


Transition 


Frequency 


Ref. 


Telescope 


HPBW a 






(MHz) 






(") 


C4H 


N — 2 — 1 /-I 5 — 5 F — 1 — 1 


1 Q054 4760 


1 


GBT 


38 




7 = 2—1 


lovj i.y . 1 jou 


2.3 


GBT 


39 


C 6 H- 


7= 10-9 


27537.1302 


4 


GBT 


26 


C 6 H- 


7= 11-10 


30290.8133 


4 


GBT 


24 


C 6 H 


7= 10.5-9.5, f,F = 11-10 


29109.6437 


5 


GBT 


25 


C 6 H 


7= 10.5-9.5, f,F = 10-9 


29109.6855 


5 


GBT 


25 


C 6 H 


7= 10.5-9.5, e,F = 11-10 


29112.7087 


5 


GBT 


25 


C 6 H 


7= 10.5-9.5, e,F = 10-9 


29112.7503 


5 


GBT 


25 


HC,N 


7 = 3-2, F = 3-3 


27292.9018 


6 


GBT 


27 


HC3N 


7 = 3-2,F = 2-l 


27294.0758 


6 


GBT 


27 


HG,N 


7 = 3-2,F = 3-2 


27294.2927 


6 


GBT 


27 


HC 3 N 


7 = 3-2, F = 4-3 


27294.3451 


6 


GBT 


27 


HC3N 


J = 3-2,F = 2-2 


27296.2334 


6 


GBT 


27 


HG,N 


7 = 4-3, F = 4-4 


36390.8861 


6 


GBT 


20 


HC 3 N 


7 = 4-3, F = 3-2 


36392.2358 


6 


GBT 


20 


HC3N 


7 = 4-3, F = 4-3 


36392.3293 


6 


GBT 


20 


HC,N 


7 = 4-3, F = 5-4 


36392.3630 


6 


GBT 


20 


HC3N 


7 = 4-3, F = 3-3 


36394.1765 


6 


GBT 


20 


HG,N 


7=10-9,^ = 9-8 


90978.9838 


6 


OSO 


42 


HC 3 N 


7= 10-9,F = 10-9 


90978.9948 


6 


OSO 


42 


HC 3 N 


7=10-9,^ = 11-10 


90979.0024 


6 


OSO 


42 


References. 


— (1) Gottlieb et al. 1 1983); (2) Gupta et al. 12007); (3) McCarthy & Thaddeus (2008); (4) McCarthy 


■ et al. (2006); (5) McCarthy & Thaddeus 


12005); (6) HC3N line predictions are based on spectroscopic data summarized by Thorwirth et al. (2000) and de Zafra 


1197 ID. 





a Telescope half-power beam-width at observed frequency. 



TABLE 3 

Molecular line measurements and anion-to-neutral ratios 



Source 


/ T M Bdv(C 6 Erf 


/ 7ksdv(C 6 H) b 


N(C 6 U-) 


N(C 6 H) 


C 6 FT/C 6 H 




(mKkms- 1 ) 


(mKkms -1 ) 


(10 10 cm" 2 ) 


(10" cm- 2 ) 


(%) 


LI 172 


6.7 (0.8) 


41.1 (1.6) 


2.4 (0.3) 


7.1 (0.3) 


3.3 (0.5) 


L1251A 


6.5 (1.0) 


43.6 (2.1) 


2.3 (0.4) 


7.6 (0.4) 


3.0 (0.6) 


L1389 


5.9 (0.8) 


27.1 (1.3) 


2.1 (0.3) 


4.7 (0.2) 


4.4 (0.8) 


L1495B 


9.6(1.0) 


141.6 (2.2) 


3.4 (0.4) 


24.6 (0.4) 


1.4 (0.2) 


L1512 


4.3 (0.4) 


26.3 (0.9) 


1.5 (0.1) 


4.6 (0.2) 


3.3 (0.4) 


TMC-1C 


13.6(1.1) 


88.1 (1.7) 


4.8 (0.4) 


15.3 (0.3) 


3.1 (0.3) 


TMC-1 CP 


41.6 (5.0) 


332.4 (9.0) 


14.7(1.8) 


57.8(1.6) 


2.5 (0.4) 


Note. — 


Uncertainties given in parentheses 


are ±68% Monte Carlo errors. 


Measurements for L1251A and L1512 supersede those of Cordiner et al. (201 1 ) 


as explained 


in Section|4] 











a Integrated intensities of average CgH 7 = 10 — 9 and 7 = 1 1 — 10 spectra. 

b Integrated intensities summed over both hyperfine components of the CgH 7 = 10.5-9.5, / transition. 
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FIG. 3. — Observed HC3N spectra (black points) with RADEX least-squares models overlaid (red curves). Frequency- switching residuals have been masked 
from the observed data. Three main hyperfine peaks are visible in the 7 = 3-2 and 7 = 4 — 3 spectra. Model parameters are given in Table |4] 
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TABLE 4 

HC 3 N RADEX FIT RESULTS 



Source 


W(HC 3 N) 
(10 13 cm- 2 ) 


v(HC 3 N) 
(km s- 1 ) 


Av(HC 3 N) 
(kms- 1 ) 


«H 2 

(10 4 cm" 3 ) 


L1172 


2.7 


2.74 


0.49 


7.5 (2.3) 


L1251A 


3.3 


-3.97 


0.31 


2.1 (0.4) 


L1389 


1.5 


-4.75 


0.32 


5.2(1.3) 


L1495B 


8.2 


7.62 


0.19 


1.1 (0.3) 


L1512 


4.2 


7.11 


0.10 


2.6 (0.9) 


TMC-1C 


6.4 


5.37 


0.18 


1.1 (0.4) 


TMC-1 CP 


19.5 


5.75 


0.32 


1.0 (0.2) 



NOTE. — Uncertainties on «h 2 (in parentheses) were calculated assuming 
±1 K uncertainty on the kinetic temperature, but may be larger due to the 
possible effects of beam dilution. 



Doppler line FWHM (Av) and Doppler velocity (v) were op- 
timized for the spe ctra of each sour ce using the MPFIT least- 
squares algorithm (Markwardt 2008). The best-fitting param- 
eters are given in Table @] The kinetic temperature was fixed 
at 10 K, then varied by ±1 K to produce the quoted error 
estimates on «h 2 - For all sources, a temperature of 10 K pro- 
duced a good fit to all five HFS components (the two weak- 
est HFS components are not shown in Figure |3). The HC3N 
Doppler velocities are shown with vertical dashed lines in Fig- 
ure [2] and closely match the central velocities of the detected 
hydrocarbons and anions. In several cases, the HC3N line 
widths are very narrow; for L1512, Av = 0.1 kms" 1 , which 
is only slightly greater than the 10 K thermal line-width of 
0.07 kms" 1 , indicating an unusual lack of turbulence, flows 
or shears in the central ~ 30" of this source. 

The C(,H~/C^H ratio is plotted as a function of «h 2 in 
Figure [4] and shows a positive correlation. There is con- 
siderable scatter a nd uncertainty on n^ , but consistent with 
previous studies (|Hirahara et al.l I1992[ [Visser et al.l I2002[ 
Launh ardt et al.ll2Q 1 01 and lKirk et al.H2.Q05l respectively), we 
find the lowest density in TMC-1, highest densities in in 
LI 172 and L1389, and an intermediate density in L1512. 
Our calculated «h, values for LI 172, L1389 and L1512 
are, however, systematically lower than previously-derived 
values by up to two orders of magnitude. This discrepancy is 
likely indicative of a problem in our method, which can be 
attributed to the effects of different degrees of beam dilution 
affecting the different HC3N line frequencies we observed. 
Given the larger 42" beam size for the HC3N J = 10-9 line 
observed with the OSO 20-m compared to the 27" GBT 
beam for the J = 3-2 line, it is theoretically possible that the 
OSO data suffers up to an additional factor of 2.6 in beam 
dilution compared with that of the GBT, which could result 
in sufficient reduction of the derived values to account 
for the observed d iscrepancies. Figure Q] (and Figure 7 of 
Pratap et al.lll997l) shows that the HC3N emission from all 
our observed sources is distributed over an area larger than 
the OSO beam, which suggests that beam dilution may not 
be so severe. However, given the lack of spatial information 
on the extent of the J = 10-9 emission on size scales less 
than the OSO beam size, it is impossible to quantify the 
importance of beam dilution in our RADEX calculations, 
so the «h, values given in Table |4] should be treated with 
caution. More accurate estimates for «h 2 could be obtained 
from HC3N J = 10-9 observations using a telescope with a 
smaller beam (such as the IRAM 30-m), which would be a 
better match to that of our GBT observations. 
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FIG. 4. — CjH anion-to-neutral ratio vs. H2 number density («Hi)- Dashed 
line shows linear least-squares fit. 

5. DISCUSSION 

5.1. C(,H~ 

Since the 2006 discovery of CeH" in TMC-1, there have 
been seven reported detections of this anion in various parts 
of the Taurus-Auriga Molecular Cloud complex (including 
our latest detections in L1495B, L1512 and TMC-1C). Thus 
far, there have been only tw o detections of interstellar anions 
outsid e of this small region (ISakai et al.ll2010t ICordiner et alj 
201 1), which highlights the importance of our new discover- 
ies of C(,H~ in the vicinity of LI 172 SMM (a prestellar core 
in Cepheus) and L1389 SMM1 (a prestellar/protostellar core 
inside Bok Globule CB17 in Camelopardalis). These new 
detections confirm that anions are indeed widespread in star- 
forming regions outside of Taurus, and that C(,H~ appears to 
be ubiquitous in dense interstellar clouds wherever its parent 
n eutral is present. 

lMverset al.1 d 19881) detected an outflow in LI 172 that 
iVisser et al.1 "(12002) identified as originating from a protostar 
(LI 172 SMM1) that lies about 30" south of our targeted po- 
sition. They also detected a dense core (presumed prestel- 
lar), about 50" SE of our target, with N(H 2 ) = 6 x 10 21 cm" 2 . 
If this column density also applies to our observed position, 
the C(,H~ abundance is 4 x 10~ 12 , which is a factor of a few 
less than in TMC-1, but an order of magnitude greater than in 
L1512, L1521F and L1544 (see ICordiner & Charnlevll20T2l 
and references therein), suggesting a relatively high absolute 
anion abundance in LI 172. 

By utilizing the available physical and chemical informa- 
tion on the interstellar clouds for which CeH" has so far 
been detected (as summarized for our observed targets in Sec - 
tion [2), it is possible to examine the C6H"/C6H ratio with 
respect to each cloud's chemical and dynamical evolution- 
ary state. The evolutionary states of the observed clouds 
are given in Table Q] and are shown by the different sym- 
bols on the plot in Figure |4] There is a clear division in 
the anion-to-neutral ratio between young, 'quiescent' and 
older, 'star-forming' cores: for the quiescent, chemically- 
young cores L1495B and TMC-1 (CP), C 6 H"/C 6 H < 3%, 
and for the star-forming cores/protostars, CeH"/C6H > 3%. 
L1495B is the most chem ically-young and quiescent of our 
sources dHirota et al.ll2004|) , with no evidence for nearby star- 
formation, and has the lowest C^H'/C^H ratio in our sam- 
ple (1.4 ±0.2%). The very low-luminosity, low-mass pro- 
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tostar (VeLLO) L1389 SMM1, on the other hand, evidently 
contains among the most evolved matter in our sample (hav- 
ing passed through the quiescent and star-forming stages), 
and has the highest anion-to-neutral ratio (CeH"/CeH = 4.4 ± 
0.8%). This value i s similar to the v alue of 4.0 ± 1.0% 
found in L1521 F bv iGupta etal.1 (120091) . which also con- 
tains a VeLLO dBourke et al.M2.Q06b . Comparing our results 
with other sources from the literature, L1527 (Sak ai et ail 
2007) contains a somewhat more-evolved Class 0/1 protostar, 
and this source has a significantly larger CeH"/CeH ratio of 
9.3 ±2. 9%. Thus, we confirm the apparent relationship iden- 
tified by Cordi ner et al.l (|201 1 ) between the evolutionary state 
of interstellar matter and its C(,H~ anion-t o-neutral rat i o. The 
quiescent cloud Lupus- 1 A observed by ISakai et all (120101) 
also fits the trend (with C 6 H"/C 6 H = 2.1 ±0.6%), as does 
the very low-luminosity protostar Chameleon MMS1 (with 
its upper limit of C 6 H _ /C 6 H < 10%; ICordiner et all 120 121) . 
The dense prestellar core LI 544 m ay deviate, howeve r, with 
a relatively low ratio of 2.5 ± 0.8% (IGupta et al.ll2009l) . 

The observed CeH"/CeH trend and its correlation with «h 2 
shown in Figure |4] can be understood in the co ntext of the 
theoretical study of Cord iner & Charnleyl (120121) . who iden- 
tified the effects of depletion on the C(,H~ anion-to-neutral 
ratio, the degree of which is related to both the density and 
chemical age of the cloud. Depletion occurs as atoms and 
molecules collide with and stick to dust grains, and thus pro- 
ceeds faster in denser media. Increased depletion affects the 
anion-to-neutral ratio in two ways: (1) the free electron abun- 
dance goes up so that more C(,H~ is produced by radiative 
electron attachment and (2) the atomic O and H abundances 
go down, thus reducing the anion destruction rate. Our CeH" 
observations agree with this theory, particularly for the quies- 
cent cores and the protostars, which are believed to lie at op- 
posite ends of the interstellar chemical/dynamic evolutionary 
path. Measured anion-to-neutral ratios, however, are not suf- 
ficiently accurate to draw a clear distinction among the (mod- 
erately evolved) prestellar cores and star-forming gas clouds 
of LI 172, L1251A, L1512 and TMC-1C. These clouds likely 
contain gas that spans a range of densities and degrees of 
chemical evolution. More accurate measurements of density, 
depletion and the CeH"/CeH ratio at high spatial resolution 
will be required in order to confirm the utility of the anion- 
to-neutral ratio as a measure of the evolutionary state of inter- 
stellar clouds. 

5.2. C 4 H~ 

Our detection of C4H" constitutes the smallest reported col- 
umn density for this molecule in any source to-date. Consis- 
tent with the trend described above for CeH", the C4H" anion- 
to-neutral ratio of (1.2 ± 0.4) x 10" 5 in TMC- 1 (CP) is about 
nine ti mes smaller than observed in LI 527 by lAgundez et al.1 
(2008). Theory regarding C4H" chemistry is less well under- 
stood than for CeH", but it seems plausible to again ascribe 
the lower value in TMC-1 (and the moderately low value of 
(8.8 ±5.3) x 10" 5 in Lupus-lA; ls"akai et al.ll2010l) . to chemi- 
cal effects resulting from the lower densities and younger evo- 
lutionary states of these sources compared with L1527. 

The very low observed C4H" anion-to-neutral ratios com- 
pared with CeH" are at variance with chemical models that 
consider the formation of C4H" to be by r adiative electron at- 
tachm ent to C4H at the rate calculated bv iHerbst & Osamural 
(2008). Recent ab-initio calculations by V. Kokoouline and 
colleagues (private communication, 2012) show that the C4H 
radiative attachment rate may be several orders of magnitude 



less than previously thought, and based on o b servat ions of 
the C4H-/C4H ratio in LI 527, lAgundez et all (120081) calcu- 
lated a C4H radiative atta chment rate a factor of 1 22 less 
than the theoretical value of lHerbst & Osamural (120081). How - 
ever, using the laboratory d ata of Eichelberger et al.l (120071) . 
ICord iner & Charnleyl (1201 2l) identified that a significant path- 
way to smaller hydrocarbon anions is via reactions of larger 
anions with atomic oxygen. Thus, an important source of 
C4H" is the reaction 

C5H- + O— >C 4 H" + CO. (1) 

As a consequence of this (and the analogous reaction 
that forms CsH~ from CgH~ + O), using the model of 
ICordiner & Charnleyl (120 12b at n H , = 10 4 cm" 3 (applicable to 
TMC-1), it is possible to set the radiative attachment rates for 
both C4H and C5H to zero and still produce an absolute C4H" 
abundance and C4H"/C4H ratio that are within a factor of a 
few of the observed values. It is therefore plausible that ra- 
diative electron attachment is not the dominant route to the 
formation of C4H" in molecular clouds. 

6. CONCLUSION 

We conducted a search for the carbon-chain anion CeH" in 
seven nearby molecular clouds selected based on their strong 
HC3N emission. CeH" and CeH were detected in all sources, 
with a mean anion-to-neutral ratio of 3.0% (standard devia- 
tion 0.92%). Combined with the four previous CeH" detec- 
tions in other sources, we deduce that anions are ubiquitous 
in the ISM wherever sufficient CeH is present to be easily de- 
tectable. 

From the combined sample of eleven known interstellar 
sources with CeH" detections , we confirm the trend identi- 
fied by Cordin er et al.1 d2011) for anion-to-neutral ratios to 
be smaller in younger, less dense environments and larger in 
older, denser environments. This i s consistent with the the- 
ory of ICordiner & CharnleyJ (120121) that anion-to-neutral ra- 
tios are sensitive to the depletion of atomic O and H, and to 
the electron density, both of which become larger in denser, 
more chemically-evolved cores. Molecular hydrogen num- 
ber densities were derived for each source based on radiative 
transfer fits to observed HC3N spectra; a positive correlation 
is observed between CeH"/CeH and «h 2 . 

We report the first detection of C4H" in TMC-1, and de- 
rive an anion-to-neutral ratio C 4 H"/C 4 H = (1 .2 ± 0.4) x 10" 5 
(=0.0012 ±0.0004%). This is the smallest anion-to-neutral 
ratio yet observed and confirms a problem with the theoreti- 
cal value for the C4H electron attachment rate, which is ap- 
parently several orders of magnitude too large. The observed 
C4H" may instead be formed by reactions of the larger hydro- 
carbon anions (CeH" and C5H") with atomic oxygen. 

Given the widespread presence of anions in space, a combi- 
nation of further dedicated laboratory, observational and the- 
oretical studies will be required in order to better understand 
their behavior and thus their full importance in astrochemistry. 
Measurements of radiative electron attachment rates for hy- 
drocarbons, and branching ratios for reactions of anions with 
atomic oxygen will be of particular importance in this regard. 

This research was supported by the NASA Exobiology Pro- 
gram and the NASA Astrobiology Institute through the God- 
dard Center for Astrobiology. We gratefully acknowledge the 
assistance of Tom Millar and Catherine Walsh during the 2010 
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